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ABSTRACT 

With the success of the Chandra and XMM missions and the maturation of gravitational lensing 
techniques, powerful constraints on the orbital structure of cluster dark matter halos are possible. I 
show that the X-ray emissivity and mass of a galaxy cluster uniquely specify the anisotropy and velocity 
dispersion profiles of its dark matter halo. I consider hydrostatic as well as cooling flow scenarios, and 
apply the formalism to the lensing cluster CL0024-I-16 and the cooling flow cluster Abell 2199. In both 
cases, the model predicts a parameter-free velocity dispersion profile that is consistent with independent 
optical redshift surveys of the clusters. 

Subject headings: Galaxies: clusters: general — hydrodynamics — cooling flows — X-rays: 
galaxies — galaxies: clusters: individual (CL0024-I-16, Abell 2199) 



L INTRODUCTION 

In the standard picture, a cluster of galaxies consists of 
a dark matter halo accompanied by galaxies and a hot, 
1-15 keV plasma. The techniques for measuring the mass 
of the halo usually focus exclusively on the motions of the 
galaxies in the cluster, on the emissivity and temperature 
of the X-ray emitting gas, or on the gravitational lensing 
of distant sources by the halo. In this paper I show that 
combining these various data can place strong constraints 
on the velocity structure of the cluster dark matter halo. 

Over the past two decades, a great deal of evidence has 
mounted that the X-ray emissivity e of the plasma for 
many clusters is well fit by the particularly simple function 
(e.g. Jones & Forman 1984; Mohr, Mathiesen, & Evrard 
1999) 



1 



-3/3 



(1) 



This function is often referred to as the /3-model in the 
literature; it has a core radius, rc, and a faint end slope, 

^- . . . 

It is instructive to review the physical basis for the form 

of the /3-model (see, e.g. Sarazin 1988). The starting point 
is the assumption that the plasma is in equilibrium with 
the dark matter potential. The general form of the equa- 
tion of hydrostatic equilibrium is 



VP = -PgV^, 



(2) 



(3) 



where P is the gas pressure, pg is the gas density, and is 
the gravitational potential. With the assumption that the 
plasma is a spherically symmetric ideal gas, the equation 
becomes 

1 d I pgkBT\ _ GM 
Pg dr \ rupp. J 

where T is the plasma temperature, M is the total mass 
inside the radius r, is the proton mass, p is the mean 
molecular weight, and is Boltzmann's constant. The 
dynamics of the dark matter halo, on the other hand, 

^In the stellar dynamics literature, which precedes the discussion 
parameter. But because the f3 is so widely used by X-ray astronomers 



are determined by the Jeans equation (e.g. 
Tremaine 1987, p. 204), 



1 d 
Pdm dr 



+ 



GM 



Binney & 



(4) 



Here 47rr^pdm = dM/dr is the dark matter density, ar 
is the dark matter radial velocity dispersion, and rj is 
the anisotropy parameter^, which describes, in an aver- 
age sense, the nature of the orbits of the dark matter 
particles. The anisotropy parameter is a dimensionless 
quantity equal to 1 — ag/a^, where is the tangential 
velocity dispersion. When rj — 1, the dark matter orbits 
are completely radial; when rj = — oo, the orbits are com- 
pletely tangential; and when = 0, the orbits are isotropic. 
Note that equation (^) is not the same as the galaxy Jeans 
equation, which involves the galaxy number density fgai in 
place of the dark matter density pdm- 

Three critical assumptions shape the derivation of the /3- 
model. The first is that the gas is isothermal, and the sec- 
ond is that the dark matter orbits are perfectly isotropic. 
Then, equating (||) and (^), 



dlnpg d\n Pdm pmpcr"^ 



din? 



dim 



knT 



(5) 



The final assumption is that the dark matter halo is nearly 
an isothermal sphere, with a density profile given by a King 
model, Pdm oc (1 + r-^/r^)^^/^. Then the above equation 
has Pg (X (1 -hrVr2)-3/5/2^ with /3 = prupa^ / (kBT). The 
isothermal emissivity, proportional to the square of the gas 
density, is then given by equation (|l|). 

There is an inconsistency in the above derivation. The 
assumption that rj = everywhere is not compatible with a 
constant velocity dispersion ar, easily verified by inserting 
a King model into equation (^), setting ry = 0, and notic- 
ing that fJr must vary with the radius. Yet the /3-model 
derived in this manner continues to provide an adequate 
empirical description of cluster emissivities. 

of X-ray clusters by many decades, it is (3 that denotes the anisotropy 
to refer to the faint-end slope of the emissivity, I must switch symbols. 
1 
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It is interesting to ask what anisotropy profiles r]{r) and 
velocity dispersion profiles o'^(r) are fully consistent with 
the /3-model taken in conjunction with auxiliary measure- 
ments, e.g. gravitational Icnsing. A further impetus for 
investigating rj is the recent development (e.g. Geller, Di- 
aferio, & Kurtz 1999) that the density profiles of some 
clusters of galaxies are well fit by a cuspy model suggested 
by Navarro, Frenk, & White (1997; NFW); see Table 0. 
Because gasdynamical N-body simulations indicate that a 
/3-model gas can coexist with an NFW dark halo (e.g. Eke 
et al. 1998), an analytic exploration of the basic equilib- 
rium properties of the anisotropy parameter is appropri- 
ate. Finally, a theoretically sound functional form for 77 is 
required for analyzing cluster velocity dispersion profiles 
using the Jeans equations (e.g., Carlberg et al. 1997). 

In this paper I develop a formalism for calculating 77 (r) 
and tT^(r) simultaneously. I assume that the cluster is 
spherically symmetric, that the gas is in hydrostatic or 
quasihydrostatic equilibrium with the dark matter, that 
dark matter dominates the cluster potential, and that the 
gas temperature is always proportional to the local dark 
matter velocity dispersion. In §0 I derive the basic equa- 
tions; in I apply the model to hydrostatic plasmas and 
the cluster CL0024-I-16; in §|[ I consider anisotropy pro- 
files within quasihydrostatic cooling flows, specifically the 
cluster Abell 2199; and in §| I summarize. 



2. DERIVATION 

The aim of this derivation is to solve the dark mat- 
ter Jeans equation (||) for the radial velocity dispersion 
ar and the orbital anisotropy parameter rj given a mass 
profile M. Because there are two unknown functions and 
only one equation, auxiliary assumptions are necessary. 
One approach is to make ry a constant (Mahdavi et al. 
1999; van der Marie 2000) and solve for ar{r, rf). However, 
N-body simulations suggest that the orbital distribution 
of the dark matter particles varies considerably through- 
out the halo (e.g. Kaufmann et al. 1999; Diaferio 2000). 
Another method involves assuming a functional form for 
the phase space distribution of the dark matter particles 
(Merritt 1985; Gerhard et al. 1998). This approach is tan- 
tamount to fixing the form of ry(r). 

Here I outline a third approach that places no require- 
ments on 77 itself. Instead I assume that the total specific 
energy of the dark matter particles is everywhere propor- 
tional to that of the plasma: T oc cr^of This scaling 
law is not arbitrary, but is grounded in observations of 
groups and clusters of galaxies. Combined X-ray and op- 
tical observations of these systems have found that the 
mean emission-weighted gas temperature scales roughly 
as the second power of the total galaxy velocity dispersion 
(Mulchaey & Zabludoff 1998; Xue & Wu 2000). These 
data show that in many clusters, the dark matter, gas, 
and the galaxies are in dynamical equilibrium, with the 
specific energy of each component proportional to that of 
the other two at all observable radii. 

Now I describe the implications of T cx ct^q^. for the 
anisotropy parameter. In a spherically symmetric dark 
matter halo, the velocity dispersion tensor is close to di- 
agonal in spherical coordinates. The local energy per unit 
mass is therefore given by the trace of the tensor, the total 



three-dimensional velocity dispersion: 



(3 - 2r,) , 



(6) 
(7) 



where ae and ct^ are the tangential and azimuthal veloc- 
ity dispersions, respectively. In spherically symmetric sys- 
tems, as = cr^ and and ry = 1 — (Xg/cr^. Then the hypothesis 
that T (X a^Q^. implies the following relationship between 
the gas temperature and the velocity dispersion: 



(8) 



where /3truo is a dimensionless constant. 

Equations (^), (^, and (||) now form a closed set, solv- 
able as follows. I begin by writing the equation of hydro- 
static equilibrium in a different form: 



dlnpgT GMnipfi 



din 'I 



knTr 



(9) 



In clusters of galaxies, the plasma radiates chiefly 
through thermal bremsstrahlung, a coUisional process 
whose total emissivity scales roughly as the square root 
of the temperature: e oc P^Vt. It follows that 

In ~ - In e — - In T -f constant. (10) 

It is then possible to reformulate equation (||) in terms of 
the observed emissivity: 



3dlnr I dine 



4 dim 



2dlni 



GMrUpfi 
ksTr 



(11) 



This relation implies that, for isothermal plasmas, there is 
a one-to-one correspondence between the measured emis- 
sivity profile and the mass profile. In other words, an 
isothermal sphere of ideal gas with a measured emissivity 
profile e(r) must have the mass profile given by setting 
dlnT/dlnr = in equation (|l): 



M(r) = 



ksTr dine 
2Gmpfi dlnr 



(12) 



In general, the plasma can have a temperature gradient: 
the profile is given by the solution of equation ([ll[): 



T(r) 



AflTTlp 



CM 



^'^dr. 



(13) 



Here the constant of integration is chosen such that T con- 
verges as r — s- oo. The temperature at infinity is 



lim 

r — *C30 



2GMmpfi /dine 



ksr 



din 'I 



(14) 



If the logarithmic derivative of the emissivity tends to 
a constant value, then as r — s- oo, (1) Too = unless 
M{r) oc r, and (2) M{r) oc r is the steepest mass profile 
allowed. A further useful result is obtained by applying 
I'Hoptial's rule to the derivative of equation |l^: 



,. dlnT dlnM , 
hm — = lim — 1. 



din I 



din I 



(15) 
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That is, T{r) cx M{r)/r asymptotically. Because the mass 
can never decrease, the temperature never falls faster than 
at large radii. 

Now I rewrite the Jeans equation for a spherically sym- 
metric system of collisionless particles: 



dlnpdmCT^ , „ GM 

— — + 2?7 = 

amr a^r 

Using equation (^), it is possible to eliminate ry: 



din pdmO-^ 
dlnr 



3 1- 



GM 



This equation has the solution 



-r pdmdr, 



'2 — ~3 



GMrpdmdr. 



(16) 

(17) 

(18) 
(19) 

(20) 



The constant of integration should be zero in most 
cases, because a nonzero constant will usually introduce a 
if'^Pdin)~^ divergence in the velocity dispersion as r — > 0. 
Such a divergence cannot be reconciled with equation (|8|) 
if the temperature is finite at r = and t] < 1. 

Two physical considerations limit the possible values of 
/3truo- The first is that cr^ > 0. Then 



172 < S/^trucCTl- 



(21) 



The second consideration is 77 < 1, equivalent to requir- 
ing that the tangential velocity dispersion ag always be 
positive. Equation (H) then implies 



2^oa (.2 kBT\ 



(J2 > 3/3tri 



(22) 
(23) 

(24) 



Any physically acceptable emissivity-potential pair must 
simultaneously meet the requirements of inequalities (21) 
and (HJ) at all radii. Only values of /?truo which meet 
these two constraints are valid. If at any two radii the two 
inequalities are in conflict, the emissivity-potential pair 
is not physically consistent with the assumptions of this 
derivation. Whenever T — > asymptotically and 2 do 
not, the above inequalities force a unique value of Ptmo- 



Ptru 



lim „ . 

1 — >oo 3(TJ 



(25) 



Finally, for a valid choice of /3truc , the anisotropy profile 
is given by the inversion of equation (||) : 



T](r) 



1 



Ptruck bT 

mppa^ 



(26) 



Thus only the shapes, and not the normalizations, of the 
gravitational potential and the emissivity profile determine 
the velocity anisotropy. 



3. HYDROSTATIC SOLUTIONS 

Here I discuss the anisotropy solutions for clusters in 
perfect hydrostatic equilibrium — i.e., assuming radiative 
cooling is unimportant. I begin with potential-emissivity 
pairs described by simple power laws. Although single 
power laws do not provide a satisfactory description of real 
clusters, they elucidate the general properties of the solu- 
tions more clearly than broken power laws, for which the 
integrals in equations (|l^)-(20) are often not expressible 
in terms of elementary functions. 

3.1. Power Law Mass and Emissivity Profiles 

Consider a cluster with mass profile M — Ar" (with 
Pdm r"""^), and a power law emissivity profile e oa r'^. 
The mass must never decrease with the radius, so a > 0. 
Then equation (O) has the form 



AGAr" 



kT{r) _ 
mpfi 3(1 — a — 1v jZ) 



(27) 



It is unphysical for the temperature to become infinite at 
large radius. Furthermore, the temperature must always 
be positive. Hence we have the constraints 



< a < 1 

3 — 3a 
V < — - — 



(28) 
(29) 



Thus, if the matter distribution behaves like a singular 
isothermal sphere (a — 1), the plasma is isothermal as 
long as its emissivity is a declining power law of arbitrary 
index. 

The radial velocity dispersion from equation (llE 



7 = 



GAr" 



2a- 
1 - a 



'1 

1 \ 7 
- 2v/3. 



- 1 



(30) 
(31) 



Again, for a singular isothermal sphere the radial velocity 
dispersion is constant regardless of the power law index of 
the emissivity. The logarithmic a = 1/2 solution is not 
permitted, because it gives ri = 3/2 at r = 0. 

Substituting equations (27) and (31) into the anisotropy 
equation (^ ) yields. 



. = ^(3 



7//3t: 



(32) 



Thus the anisotropy parameter 77 is always constant for the 
permitted power law mass and emissivity profiles. Note 
that the dark matter orbits are perfectly isotropic (?/ = 0) 
only if the ratio of the dark matter velocity dispersion to 
the gas temperature is 



/^truc — 



37 



16- 8a 



(33) 



Not all values of /3truc are allowed. The requirements 
that cr^ > and 77 < 1 place the following bounds on 

Ptrue ' 



^ < Aruc < 7 if a < 1/2; 
— 8a 4 



< /3truo < 



7 



8 -8a 



if a > 1/2. 



(34) 
(35) 
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Fig. 1. — Limits on /3true, the ratio of the total dark matter velocity dispersion to the gas temperature, for mass and emissivity profiles 
M oc r" and e oc r". Two different views of the same plot are shown for clarity. The meshed surfaces indicate the upper and lower limits on 
Ptrue- The smooth surface shows the required value of /3true such that the velocity anisotropy 77 = 0. 
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Figure |l| depicts these boundaries, as well as the r/ = 
criterion. For a < 1, /3truc always has firm upper and 
lower limits, but for the singular isothermal sphere, arbi- 
trary large values of /Struc are allowed. Within the upper 
and lower limits, the anisotropy rj takes on the full range 
of allowed values, {— oo, 1). 

A physical explanation of overall tilt of the planes in 
Figure lies in the slope of the emissivity profile. When 
the gas distribution is relatively extended {i/ ^ —1), the 
gas is on the average hotter, and the dark matter to gas 
energy ratio /3true is smaller, than when the gas is rela- 
tively concentrated {ly < —1). Thus the mean permitted 
value for Ptmc increases as the emissivity profile steepens. 

The asymmetry in Figure |l| between the boundaries for 
a < 1/2 and a > 1/2 is related to the differences in the 
circular velocity profile of the potential, Vc = \J GMjr. 
In systems with a. > 1/2, Vc is nearly constant, and am- 
ple kinetic energy is available at all radii to particles with 
tangential orbits. It is therefore statistically difficult to 
underpopulate the tangential orbits and achieve 77 w 1. 
Thus /3truo must be inordinately large to achieve a signif- 
icant radial anisotropy: a large amount of kinetic energy 
must reside in the radial orbits for them to dominate the 
total velocity dispersion. 



different from the mass profile of the King sphere (Table 
|l]). The reason for this discrepancy is that the classical 
derivation of the /?-model emissivity assumes that the King 
sphere is isothermal (see the Introduction). But the King 
sphere is only an approximation to the true nonsingular 
isothermal sphere (Binney & Tremaine 1987, p. 228). For 
a /3-model gas to be truly isothermal, it must be embedded 
in the mass profile given by equation (36). 

It is then straightforward to derive the orbital proper- 
ties of the spherically symmetric dark matter halo using 
the equations derived in §S. They yield Ptiue = f3, and 



3l3kT 



3(1 a;2)2 tan"^ x 
x{3 + x^) 



(37) 



The central velocity dispersion is cr^(O) = (3kT/{mp^), 
and as r — > 00, cr^ ~ 1/r. Because the radial velocity dis- 
persion eventually vanishes, the orbits at r = cx3 must be 
completely tangential to keep the total velocity dispersion 
constant, and therefore 77 — > — cx). The required anisotropy 
parameter is independent of f3: 



3 2x^(3 + x^) 

^ 2 ^ i{l + x^fia.n-^x-x[3 + x'^)' 



(38) 
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Fig. 2. — Anisotropy profile (top) and radial velocity dispersion 
profile (bottom) for the isothermal /3-model. 

3.2. Isothermal [3 -models 

The mass and emissivity profiles of real clusters are bet- 
ter described by broken power laws. The /9-model emissiv- 
ity, which has e oc (1 -I- x^)"'^'', where x = r/rc, fits many 
clusters without cooling flows, and is commonly taken to 
describe an isothermal gas. If the temperature is constant, 
and the plasma is in hydrostatic equilibrium with the grav- 
itational potential, then according to equation (|l2|) the 
mass profile has to be 



M(r) 



3/3fcrrc x^ 
mp^G 1 -I- 



(36) 



Figure y shows the velocity dispersion and anisotropy pro- 
files for the isothermal /3-model. 

3.3. Non-isothermal [3-models 

Although many clusters of galaxies contain an isother- 
mal plasma, a perhaps larger fraction exhibit temperature 
gradients. For example, all but 3 of the 24 clusters ex- 
amined with ASCA by Markevitch et al. (1998) have pro- 
jected temperature profiles which decline with radius. 

A temperature gradient results whenever a /3-model 
emissivity is in equilibrium with a mass distribution dif- 
ferent from equation (|3^). Here I calculate the anisotropy 
solutions for the /3-model together with the various profiles 
listed in Table |. 

Figure |^ shows the equilibrium temperature profiles 
given by equation ( |l3| ) for (i = 0.65 and (3=1. Here 
the temperature always rises towards the core, but only 
because is the same for both the emissivity and the mat- 
ter density. Were Tc to be different for the two profiles (as 
in §3.4 below), the temperature could still fall towards the 
center. Note that the King approximation to the isother- 
mal sphere, from which the /3-model is derived, has con- 
stant temperature only within several core radii; at larger 
radii the plasma temperature drops steeply. Conversely, 
the dark matter potential of a singular isothermal sphere 
causes the /3-model gas to have a temperature gradient in- 
side the core, but a constant temperature at large radii. 
As one would expect, the gas in a cluster with (3 = 0.65 is 
on the average hotter than the gas in a cluster of the same 
mass but with (3=1. 

The next step is to decide which values of (3trnc, the ra- 
tio of the total dark matter energy to the gas energy, are 
appropriate. Figure ^ shows the upper and lower limits 
for the singular isothermal sphere from equations (^ij) and 
([2^). Because no single value of Ptme can simultaneously 
satisfy the limits at all radii, the SIS and the /3-model are 
an incompatible potential-emissivity pair. 
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Table 1 
Mass Profiles Considered 



Profile p{x) M{x) 



SIS ^ 

King (l + a;2)~3/2 \n {x + ^/YTll^) ~ x / ^/YTli? 

NEW x-i(l + a;)-2 In (1 + x) - a;/(l + x) 

Hernquist x^^(l + a;)^'^ a;^(l + a;)^^ 

Plummer (1 + a;^)"''^/^ + a;2)-3/2 



Note. — The density and mass profiles are given without normalization and in terms of the scaled radius x = r/r^. 
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Fig. 3. — The equilibrium temperature profile for two /3-model 
emissivity profiles and various gravitational potentials. The mass 
profiles are normalized so that they all contain the same amount of 
mass within r = rc- 

However, as Figure ^ shows, the rest of the mass pro- 
files in Table (El) have admissible solutions when combined 
with a /S-modelgas. In contrast to the scale- free solutions 
in § p.l| , equations (|l|) and (|2|) allow one unique value of 
/3true for each potential: 



lim -\. 



(39) 



In general, the function /3truc(/3), which matches the ob- 
served emissivity, characterized by /3, with the physical 
energy ratio Aruo, is not analytic. However, the integrals 
in equation (n9|) and (pO) do have analytic representations 
for a /3-modelgas erriEedded in a Plummer sphere. The 
exact relations are 



T{x) oc 



3(4/3+ l)v^rT^ 



(40) 



/^true — 
(7r{x) OC 



4/3 -fl 



6\/r 



77(r) = 0, 



(41) 

(42) 
(43) 



where x = r /r^. 

The Plummer model is remarkably well matched to the 
/3-model gas emissivity. Its velocity dispersion is always in- 
dependent of /3, and the dark matter orbits are completely 
isotropic. This result is surprising and counterintuitive, 
because the /3-model emissivity is historically derived from 
the isotropic King model. But as I have discussed above, 
that derivation assumes that the perfectly isotropic King 
sphere is isothermal, an assumption which breaks down 
well outside the core. In contrast, though it is not isother- 
mal, a /3-model gas embedded in an 77 = Plummer sphere 
satisfies the equations of hydrostatic equilibrium perfectly. 

For the King, NFW, and Hernquist models Arue(/3) 
is also roughly linear, with the general property that 
/3truo < /3 in each case. The resulting radial velocity dis- 
persion profiles, given by equation (Eq), and anisotropy 
profiles, given by equation (]3^), appear in Figures || and 
|7| respectively. For the King profiles, (3 correlates strongly 
with the orbital anisotropy at infinity; (3—1 gives almost 
completely radial orbits at infinity, (3 = 0.65 results in or- 
bits with a greater tangential component, and /3 = 3/4 
gives perfect isotropy at all radii. 
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Fig. 4. — Upper and lower limits for /3truoi the ratio of the dark 
matter energy to the gas energy, for a /3-model gas embedded in a 
singular isothermal sphere. The figure shows that a constant /3truc 
is not allowed. 
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Fig. 5. — Upper and lower limits for /3truc, the ratio of the dark matter energy to the gas energy, for a /3-model emissivity combined with 
the various mass profiles listed in Table hi. Solid lines show /3 = 1, and dashed lines show 13 = 0.65. 
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Fig. 6. — Radial velocity dispersion profiles for the various potentials in Table y when they contain a /3-model gas. Solid lines show /3 = 1, 
and dashed lines show /3 = 0.65. 
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The anisotropy profiles of tlie NFW and Hernquist mod- 
els show a particularly interesting structure. Both profiles 
have a small radial anisotropy at r = 0, attain a peak of 
77 ~ 0.75 near r — r^, and then decline towards a final, 
positive value. 

The King and Plummer spheres exhibit featureless 77(7"), 
whereas the NFW and Hernquist anisotropies have a max- 
imum value. The physics behind this difference lies in the 
behavior of the dark matter density as r — > 0. The NFW 
and Hernquist profiles have a cuspy pdm oc for r ^ Tc, 
whereas the King and Plummer spheres have finite den- 
sity at r = 0. If the gas is in perfect hydrostatic equilib- 
rium with the cuspy models, its temperature profile rises 
as J — > 0. For the dark matter to gas energy ratio /3true to 

-^2(3-27?) 



be constant, the total velocity dispersion a^^^ 



must then rise towards the center as well. But erf is also 



coupled to the tem 
Jeans equation (ll(; 
u\ from equation 




iperature and mass distribution via the 
For r < Tc, the temperature term 
dominates. Because of the density 
cusp dX r < Tc, (t{ itself falls more slowly than the tem- 
perature, and an increasing radial anisotropy is required to 
compensate, keeping /3truc constant. At larger radii r > rc, 
however, the subtractive, gravitational term ct| competes 
with the temperature term and causes cr^ to fall rapidly. 
A larger tangential velocity component is then required to 
maintain a constant dark matter to gas energy ratio, and 
r]{r) falls. 

Remarkably, the NFW and Hernquist anisotropies in 
Figure ^ are similar to the dark matter anisotropy profiles 
in N-body simulations without gas (Cole & Lacey 1996; 
Kauffman et al. 1999; Diaferio 1999). In these simulations, 
r}(r) is positive at r = 0, achieves a maximum near the 
virial radius, and then decreases to negative values. The 
close correspondence between the results of this derivation, 
which rely closely on the structure of the X-ray emitting 
gas in the dark matter halos, and the N-body simulations, 
which do not include gas physics, is encouraging. 

3.4. Application to CL0024+16 

Now I use X-ray imaging, gravitational lensing, and 
an optical redshift survey to test the simple hydrostatic 
formalism developed above. I derive unique velocity dis- 
persion and anisotropy profiles by combining X-ray and 
lensing data for CL0024-f 16. I calculate the shape and 
normalization of the dark matter velocity dispersion pro- 
file without recourse to any free parameters. This profile, 
once projected, can be directly compared to the indepen- 
dent galaxy velocity dispersion profile, providing a power- 
ful consistency test of the complete dynamical model. 

CL0024-I-16, a z = 0.39 cluster, has ROSAT High Res- 
olution Imager X-ray data (Bohringer et al. 2000), an 
ASCA temperature measurement (Soucailet al. 2000), and 
a mass profile from strong lensing (Tyson, Kochanski, & 
Dell' Antonio 1998). These observations independently fix 
the X-ray emissivity and the mass profile. Furthermore, a 
catalog of optical redshifts for 138 galaxies in the field of 
the cluster is available (Dressier et al. 1999). 

The Tyson et al. (1998) surface mass density is charac- 
terized by a soft core within the lensing radius n: 



T.{R) = Kr 



(l+7i?Vr-?) 



(44) 



-(H-i?2/r-f)2-7' 
where R is the projected distance from the cluster center, 



Ki = 79OO/iioo-^^0 pc^^ is the central surface density, the 
slope 7 = 0.57, and r; = 35/iioo~^ kpc, with the Hubble 
Constant ffo = lOO/iioo km s~^ Mpc~^. Outside the lens- 
ing radius r; the mass profile resembles an NFW profile 



with a core radius rc — 311/iioo ^ kpc: 



S(i?) = 



M 



sec ^R/Tc 
2^r?(i?2/r?-l) r ^R2lrl - 1 



1 - 



(45) 



Here M is determined by requiring continuity with equa- 
tion (^). Thus the density profile of CL0024-I-16 resem- 
bles that of a NFW halo, but lacks a central singularity. 

The deprojected density profile is given by (Binney & 
Tremaine 1987, p. 205): 



Pdm{r) 



1 



dR 



TT Jr dR Vi?2 - r2 

This integral is well approximated by the function. 



(1 + r/ro)^ 



X IO^^Mq Mpc 



-3 



(46) 



(47) 



where tq = 0.0725ft.ioo~^ Mpc. The resultant projected 
mass profile is better than 12% accurate everywhere within 
the lensing arcs, commensurate with the « 15% error in 
the Tyson et al. (1998) parameters. The X-ray emissivity, 
on the other hand, is well described by a /3-model with 
core radius 33/iioo~^ kpc and (3 — 0.475 (Bohringer et al. 
2000). 

Because the central cooling time of the CL0024-f 16 is 
approximately the same as the age of the cluster, it should 
not contain a cooling flow (Soucail et al. 2000). Therefore 
the hydrostatic analysis discussed in the previous sections 
is appropriate. I take the mean molecular weight to be 
^ = 0.6. Note that of all the properties derived by the 
model, only the temperature depends on the choice of /i. 
Figure || shows the results. 

According to the model, the central temperature of the 
cluster should be ~ 9.2 keV. The emission-weighted tem- 
perature within the ASCA aperture should be 12 keV; 
the Soucail et al. (2000) spectroscopic analysis gives a 
90% confidence interval of 3.6-10.6 keV. Thus the pre- 
dicted plasma temperature is consistent with but some- 
what above the measurement. The discrepancy is not 
tied to the present analysis, but is an intrinsic prop- 
erty of the lensing and the X-ray data, which give only 
marginally consistent mass measurements for the cluster 
(Soucail et al. 2000). The lensing mass measurement might 
be an overestimate because of contamination by large scale 
structure projected along the line of sight (Metzler et al. 
1999). The X-ray measurement is also suspect, because 
(3 < 0.5 in theory yields an infinite X-ray luminosity for 
the cluster. 

Given the X-ray and lensing data, however, the calcu- 
lations show that the dark matter to gas energy ratio is 
/3truo = 0.561. Then it is straightforward to derive the 
anisotropy profile 77 (r) and the radial velocity dispersion 
profile (T^ir) through equations (p9|)~(p0|). The line-of- 
sight dark matter velocity dispersion profile is (Binney & 
Tremaine 1987, p. 208) 



i?2 

1 - 77 — 



i?2 



(48) 



11 




0.5 



0.1 1 
r (h-Jo Mpc) 




-0.5 - 



0.1 1 
r (h-Jo Mpc) 



7 1000 



800 




2500 
2000 

1500 



1000 



500 



n — I I I I I I 1 1 — I — I I I I I 



XVN = 1.5 
q = 0.15 ^ 
No free parameters 

(d) 



J 



J I I 



0.] 1 
r (h-Jo Mpc) 



0.1 1 
R (h^Jo Mpc) 



Fig. 8. — Derived properties of CL0024+16, and a comparison with the optical spectroscopic data: (a) predicted plasma temperature (solid 
line) and mean integrated emission-weighted temperature as observable by ASCA (dashed line); (b) predicted orbital anisotropy profile; (c) 
predicted three-dimensional radial velocity dispersion; (d) predicted line-of-sight velocity dispersion profile along with the measured galaxy 
velocity dispersion (Dressier et al. 1999). 
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The theoretical ap{r) appears in Figure || alongside the 
line-of-sight velocity dispersion profile of the galaxies. To 
calculate ap{r) from the Dressier et al. (1999) data, it is 
first necessary to determine the well-mixed cluster mem- 
bership. I first select all Dressier et al. (1999) galaxies 
in the CL0024-H16 field with redshifts 0.381 < z < 0.404, 
corresponding to a ±2500 km s~^ line-of-sight velocity rel- 
ative to the cluster rest frame. Then I apply the DEDICA 
algorithm (Pisani 1993) to determine the cluster mem- 
bership. DEDICA uses Gaussian kernels to arrive at a 
maximum-likelihood estimate of cluster membership, in 
this case 104 galaxies. Finally, I measure ap{r) by gath- 
ering the members into 10 bins of 10-11 galaxies, and cal- 
culating cTp and its associated errors using standard boot- 
strap analysis. 

The simple model discussed above well describes the ve- 
locity dispersion profile of CL0024-I-16. The statistic 
is 15 for 10 degrees of freedom, yielding an acceptable fit 
quality g = 0.15 (Press et al. 1993, p. 660). The agreement 
is remarkable given that no free parameters are allowed in 
the comparison. 

4. MULTIPHASE COOLING FLOW SOLUTIONS 
4.1. General Properties 
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Fig. 9. — The dark matter to gas energy ratio, /3truo, for cooling 
flow models with A: = 1, 2, 3 and faint-end emissivity slopes U2. 

So far this discussion has focused on a intracluster 
plasma in hydrostatic equilibrium with the gravitational 
potential. This scenario provides a good description of 
most clusters, at least for gas removed more than « 100- 
200 kpc from the cluster center. However, the plasma in 
many clusters has a cooling time significantly shorter than 
the Hubble time within r w 100 kpc, an indication that an 
inward flow may be present. It is worth exploring how the 
dark matter anisotropy changes within the cooling regime. 

Here I apply the formalism derived in §H to clusters 
which contain a cooling flow. In this case, hydrostatic 
equilibrium no longer applies, and equation (ph must be 



replaced with a set of cooling flow equations. I adopt the 
multiphase cooling flow equations discussed by Thomas 
(1998), who examines a quasihydrostatic, spherically sym- 
metric flow with an emulsion of comoving but thermally 
isolated density phases. The Thomas (1998) scenario, 
summarized briefly in the Appendix, is a refinement of 
the classical multiphase models of Nulsen (1986). 

The cooling flow equations have a free nondimensional 
parameter, fc, which together with the gas emissivity e 
completely specifies the shapes of the flow temperature, 
the mass deposition rate, and the gravitational potential. 
In the Thomas (1998) model, k describes, in an average 
sense, the nature of the phase mixture in the flow, and can 
have values from 1 to oo. Emulsions with k = oo include 
phases of arbitrarily low density; those with fc ~ 1 allow 
only a small range of densities at each radius. Thomas 
(1998) argues that the fc ~ 1 models are preferable because 
they reproduce observations of many cooling flow clusters 
where the emission- weighted temperature drops as r — > 0; 
on the other hand, the k — oo solutions have increasing 
T(r) as r — > 0, in contrast with the observations. 

Once k and the emissivity e are specified, the resulting 
temperature and mass profiles constrain the anisotropy 
profile of the dark matter halo through equations (18)- 
(|2^). I will not use the /3-model emissivity for the cool- 
ing flow calculations, first because it does not accurately 
match the observed emissivity profiles of cooling flow clus- 
ters, which are often cuspy as r — > 0, and also because it 
is inconsistent with all but the k ~ oo flows. Instead I use 
the broken power law 



e{x) cx {x"' -hx"^)" 



(49) 



where x = r/rc and 1^2 > vi- An emissivity profile 
with vi = resembles a traditional /3-model profile with 
p = ^2/6. See the Appendix for a detailed calculation. 

Again, the first step in deriving the anisotropy profile 
is calculating /3true, the constant dark matter to gas en- 
ergy ratio. Figure || shows the results. Here Ptmo depends 
strongly on the faint-end slope of the emissivity, 1^2 ■ When 
U2 ^ 3, Ptrue is completely independent of vi; for 1^2 ^ 3, 
a slight perturbation exists as a function of vi. 

One interesting property of the solutions is that /3truo 
never exceeds 0.65, because the multiphase equations do 
not allow larger values. The fc ~ 1 models, which best 
match real clusters (Thomas 1998), allow the largest range 
of /3truc; the k = 00 models always have 0.4 < Ptme < 0.5. 

Figure 10 shows the resulting anisotropy profiles. A 
wide variety of profiles occur depending on the specific 
value of 1^1, z/2, and k. However, all r]{r) begin positive 
and decline. When the emissivity is shallow — 1^2 = 2, for 
example — the anisotropy tends towards negative values at 
larger radii. However, for reasonably steep emissivity pro- 
files, the anisotropics stay constant or increase at large 
radii, and exhibit a substantial radial bias outside the core. 

Although the cooling flow solutions extend to r = cx), 
they have the greatest validity inside the cooling radius. If 
the measured emissivity of a cluster drops off very steeply 
at large radii (1/2 ^ 4) , no multiphase cooling flow solution 
can properly account for the gas physics at all radii. In 
such cases it is best to regard the true anisotropy parame- 
ter as resembling the cooling flo w so lution within the core, 
and the hydrostatic solutions (§3.3) for r > re- 
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Fig. 10. — Anisotropy profiles for dark matter lialos containing quasifiydrostatic cooling flows with the indicated /c-parameters. The 
emissivity profile characterized by ui and U2 is given in equation (f49). Also shown is /3truc(fc, i^i, 1^2)- 
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11. — Surface brightness plus background for the ROSAT 
PSPC observation of Abell 2199. The error bars show the uncer- 
tainty due to photon noise plus a 5% additional error to account for 
uncertainty in the location of the bins. The solid line shown is the 
best fit combination of a constant backsmund plus the projection 
of the broken power law from equation (Ud) along the line of sight, 
with ui = 1.35, V2 = 3.41, and rc = 0.059hioo~^ Mpc. 



4.2. Application to Abell 2199 

Here I test the cooling flow anisotropy model by apply- 
ing it to X-ray and optical observations of Abell 2199. I 
predict the shape and normalization of the projected ve- 
locity dispersion profile from the X-ray data, and compare 
it directly with independent optical spectroscopy. Here the 
cooling flow model determines the shapes of the mass and 
temperature profiles from the X-ray emissivity, and no ad- 
ditional lensing mass is required. An absolute temperature 
measurement provides the normalization. 

Abell 2199 is a z = 0.03 cluster that contains a mod- 
erate cooling flow, with a mass deposition rate M k, 
I5OM0 yr^^ (Peres et al. 1998). The Einstein observa- 
tory emission- weighted temperature is 4. 5 ±0.2 keV within 
2^-100^^ Mpc (David et al. 1993). Velocity measurements 
for 98 galaxies in the field of the cluster are also available 
(Hill & Oegerle 1998). 

To determine the shape of the X-ray emissivity, I exam- 
ine a publicly available ROSAT PSPC observation of Abell 
2199, with sequence identification RP800644N00, and a 
41000 s exposure time. I use the MIDAS/EXSAS soft- 
ware package (Zimmermann et al. 1993) to calculate the 
surface brightness profile within 0.5/iioo^^ Mpc. I then fit 
a constant background plus the surface brightness profile 
given by the projection of equation (49) along the line of 
sight. A deprojected emissivity with vi = 1.35, 1^2 = 3.41, 
and 



= 0.06/i]^QQ Mpc provides an excellent fit with a 
reduced chi squared equal to 0.6 (Figure P). The errors 
on the fitted parameters are less than 10%rNote that the 



innermost bin (not shown) is omitted from the analysis to 
minimize smearing of the emission due to the PSPC point 
spread function. 

Once vi and V2 are specified, the shapes of the tem- 
perature and mass profiles are fixed by the cooling flow 
equations. The David et al. (1993) temperature measure- 
ment fixes their normalization, and it is possible to fully 
calculate the velocity dispersion and anisotropy proflles. 
The line-of-sight velocity dispersion profile is then given 
by equation (Eq). To measure the profile independently 
from the Hill & Oegerle (1998) data, I apply the same 
procedure as that described in § |3.4| . The results of the 
analysis for both the fc = 1 and k — 2 solutions are shown 
in Figure (H). For fc = 1 Ptmc = 0.62, and for fc = 2 
Aruo = 0.57. The profiles are steeper for larger values of 
fc as r ^ 00. 

Once again, this simple formalism provides an ade- 
quate description of the measured velocity dispersion pro- 
file without recourse to any free parameters. The reduced 
chi squared is 2.2, providing an acceptable fit for 9 degrees 
of freedom. Solutions with fc > 1 are excluded by the data, 
because the velocity dispersion profile falls too quickly, af- 
firming the intuition of Thomas (1998) that multiphase 
cooling flows with a broad density phase distribution are 
less likely to be consistent with real clusters. 

5. CONCLUSION 

I use the X-ray emission from hot gas embedded in a 
dark matter halo to break the degeneracy between the 
anisotropy parameter 77 and the system mass. Given an 
assumed mass profile and and X-ray emissivity, the model 
simultaneously determines the anisotropy profile and the 
radial velocity dispersion profile. The resultant ri{r) for 
the hydrostatic case are consistent with anisotropy profiles 
from N-body simulations without gas. With the addition 
of central cooling fiows, the formalism predicts a rich va- 
riety of anisotropy profiles within the cooling radius, most 
of which, however, contain a substantial radial anisotropy 
at the cluster center and at large radii. 

I apply the formalism to two sets of observations of clus- 
ters of galaxies: the lensing cluster CL0024-I-16 and the 
cooling fiow cluster Abell 2199. In both cases, the veloc- 
ity dispersion profiles predicted by the anisotropy models 
are consistent with independently measured optical data. 
The agreement is particularly encouraging because no free 
parameters are allowed in the comparisons. 

The warning of Binney & Tremaine (1987, p. 208), that 
for spherical systems, "radically different models can be 
consistent with both the Jeans equations and the obser- 
vations," is particularly applicable here and should not be 
neglected. This paper shows, however, that a combination 
of the X-ray emissivity and other auxiliary data, including 
the lensing mass or the gas temperature, provides a use- 
ful and self-consistent starting point for constraining the 
orbital structure of dark matter halos with gas. 

This research was supported by the Smithsonian Insti- 
tution. I thank Margaret Geller, Joseph Mohr, and the 
anonymous referee for suggestions that improved the paper 
considerably. I am grateful to Lars Hernquist, Kathleen 
Kang, and Daniel Koranyi for useful discussions. 
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Fig. 12. — Derived properties of Abell 2199 for cooling flow solutions with fe = 1 (solid lines) and k = 2 (dashed lines): (a) predicted plasma 
temperature; (b) predicted orbital anisotropy profile; (c) predicted three-dimensional radial velocity dispersion; (d) predicted line-of-sight 
velocity dispersion profile along with the measured galaxy velocity dispersion (Hill & Oegerle 1998). 
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APPENDIX 



DERIVATION OF THE ANISOTROPY PROFILE FROM THE COOLING FLOW EQUATIONS 

Thomas (1998) derives equations that describe a quasihydrostatic, spherically symmetric flow consisting of an emulsion 
of comoving but thermally isolated density phases. The relevant variables are the mass, M, the temperature, T, the 
emissivity, e, the mass accretion rate, M, and the following dimensionless quantities: 



_ GMimip 
" 2rkBT ' 
1 dlnM 



k d\nr 
dlnM 



d\ni 



(Al) 

(A2) 
(A3) 



where ks is Boltzmann's constant, and fc is a dimensionless parameter that characterizes the range of densities present at 
each radius. The A: = 1 models possess a minimum density at each radius and are the least extended. The k = oo models 
include phases of arbitrary low density, and are the most extended convectively stable distributions. 
With the assumption that the cooling function is dominated by thermal bremsstrahlung (A cx pz 

and that the 

ratio of specific heats 7 — 5/3, the Thomas (1998) steady-state equations have the solutions 



^3,5/7 



(20/21 + fc)/r2e5/7rfr' 
dine" 



dim 



5 3 

Rearranging the equations to solve for the temperature yields 



dlnE 
dlnr 



(A4) 
(A5) 
(A6) 



T cx e^/^ 



20 , 

21+^ 



rh'^'dr 



(AT) 



where the expression in square brackets approaches 1 as A; —> 00. It is therefore evident that for all reasonable emissivity 
profiles, the k — 00 solutions have dT/dr < everywhere. 

Once the temperature is fixed, it is possible to constrain /3true, the ratio of the total dark matter kinetic energy to the 
gas energy. Because A-Kr^p^^ = Mx, the velocity dispersion equations (pj|)-(pO|) become 



ks 



mpiirTT,x 



4 



ks 



3/3trucT^Sxdr-, 
2T^j:^xdr. 



(A8) 
(A9) 



In cases where 3> kBT / (mpfi) at large radii, /3truo once again takes on a unique value: 



/Jtruo = lim 



If the integrals in equations ([Aq)-( A9) diverge as r ^ cx3, then, applying I'Hopital's rule, 
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3 r—*oo 
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^21 



20/21 +fc 
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dine 



d\nr V20/3 + 7fc 
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(AlO) 

(All) 
(A12) 



If, on the other hand, the integrals converge, then Ptmc must be determined numerically. After this, direct application of 
equations (ttq) and (Gq) will yield the radial velocity dispersion profile ctj. and the anisotropy profile r]{r). 



